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Abstract — A suitable methodology is needed to determine the 
optimallocations of PMUs in a power system. In addition to its 
abilityto measure voltage and current phasors, a state-of-the-art 
PMU may include other features such as protective actions. The 
scopeof the present paper is limited to the use of PMUs for state 
estimation. 

A power system is called completely observable onlywhen all 
of its states can be uniquely determine. 

The Phasor Measurement Unit (PMU) is considered to be 
one of the most important measuring devices in future for the 
power system. The distinction comes from its unique ability to 
provide synchronized phasor measurements of voltages and 
currents from widely dispersed locations in an electric power 
grid. The commercialization of the global positioning system 
(GPS) with accuracy of timing pulses in the order of one 
microsecond makes possible the commercial production of 
PMUs. 

Traditional state estimators were non-linear in nature takes 
very high computation times. Hybrid state estimator with a 
measurement model combining the state-estimate from the 
classical state estimator with the direct state measurement from 
the PMUs was linear and the solution is direct and non-iterative. 
In Hybrid state estimator by addition of even a few PMUs can 
significantly increase the accuracy of the system state. Previous 
PMU placement research has focused primarily on network 
topology, with the goal of finding configurations that achieve 
full network observability with a minimum number of PMUs 
due to high cost of PMUs. The cost associated with these devices 
has fallen over the last several decades as computer & GPS 
technology has improved while the cost of PMUs has decreased. 

Index Terms —PMU, GPS, Measurement Unit, Electric 
Power 


I. Introduction 

Rising dependence of human society on electrical energy is 
associated with a desire for reliability of the highest order in 
power system. Over the years, we have grown to depend on 
the near perfect reliability of these systems that have become a 
necessary part of our everyday lives. All of our household 
appliances, communication devices, and almost all of our 
tools ranging from construction sites to our offices require 
electricity for operation. It is not as if we assume electricity 
will always be available, it is that we believe electricity will 
always be available. 

State estimation is facet of electric power engineering that has 
evolved out of these needs. State estimator uses information 
from various measurement units and monitoring systems to 
estimate the states of a power system. Until recently, it was 
not possible to measure phase angles of the bus voltages in 
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real time due to the technical difficulties in synchronizing 
measurements from distant locations. The advent of phasor 
measurement units (PMUs) alleviated this problem by 
synchronizing the voltage and current waveforms at widely 
dispersed locations with respect to a Global Positioning 
System (GPS) clock [1]. 

A suitable methodology is needed to determine the optimal 
locations of PMUs in a power system. In addition to its ability 
to measure voltage and current phasors, a state-of-the-art 
PMU may include other features such as protective actions. 
The scope of the present paper is limited to the use of PMUs 
for state estimation. A power system is called completely 
observable only when all of its states can be uniquely 
determined [2], [3]. 

The major objective of this paper is to outline what is being 
done in the field of state estimation using PMU data 
exclusively by presenting the work completed. The end goal 
of the project is to development of the linear state estimation 
algorithms for a three phase tracking linear state estimator and 
then use this state estimator to drive several other applications 
control center including instrument transformer calibration 
and intelligent islanding. The topology processor, which will 
take breaker status and line current phasors and determine the 
most up-to-date topology of the network, and the actual three 
phase linear state estimator which will use PMU data to 
determine the state of the system . 

This paper discusses the background information on the topic 
of three phase linear state estimation and then presents the 
work done towards developing the state-estimator. It begins 
with a brief history of the significance of power system state 
estimation followed by the presentation of the traditional form 
of state estimation. It continues with the development of the 
linear state estimation equations and then investigates the 
planned to design mathematical algorithms for a three phase 
linear state estimation and the associated research and 
software development. 

II. PHASOR MEASUREMENT UNIT 

Phasor Measurement Unit (PMU) - a device which by 
employing widely used satellite technologyoffers new 
opportunities in power system monitoring,protection, analysis 
and control. PMUs facilitate innovative solutions to 
traditional utility problems andoffer power system engineers a 
whole range of potentialbenefits, including precise estimates 
of the power systemstate can be obtained at frequent intervals, 
enablingdynamic phenomena to be observed from a 
centrallocation and appropriate control actions taken. 
Postdisturbanceanalysis will be much improved with 
theprecise snapshots of the system states through 
GPS synchronization. Advanced protection based upon 
synchronized phasor measurements could be implemented, 
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with options for improving overall systemresponse to 
catastrophic events. Advanced control usingremote feedback 
becomes possible, thereby improvingcontroller performance. 
Bock diagram of a PMU shows in figure 1. In the following 
figure shows anti-aliasing filters is meant to filter out from the 
analog input waveform frequencies more than the Nyquist 
rate. The other block is phaselocked oscillator (PLL) converts 
the GPS (Global Positioning System) pulse per second into 
the sequence of high-speed timing pulses used in waveform 
sampling. The microprocessor performs the FFT Phasor 
calculations. 



Figure 1: Block Diagram of a PMU 


PMU measures voltage and current phasor in a power system. 
Synchronism among phasor measurements is achieved by 
same time sampling ofvoltage and current waveforms using a 
commonsynchronizing signal from the global positioning 
satellite [ 8 ] - [ 10]. PMUs have higher accuracy than 
conventionalmeasurements. They reduce effects of time-skew 
amongmeasurements, useful for many other applications such 
as system protection, control and stability assessment, 
aidtopology error identification, parameter error detectionand 
correction and improve accuracy of state estimation. 

The Optimal placement of PMU becomes animportant 
problem to be solved in power system stateestimation. The 
PMU placement problem is formulated asa binary integer 
linear programming, in which the binarydecision variables (0, 
1) determine whether to install aPMU at each bus, while 
preserving the system observability and lowest system 
metering economy. It isneither economical nor necessary to 
install a PMU atevery node of a wide-area interconnected 
network. 

The cost of a PMU depends on a number offactors, including 
the number of measuring terminals(channels), CT and PT 
connections, power connection,station ground connection, 
and GPS antenna connection. 

The main purpose of performing PMU placementproblem is 
to minimize the number of installed PMUs, so that for an 
n-bus system the optimization problem isgiven as: 



Subject to /(X) > 1 

Where the installation is cost of the PMU at bus i. 


Assume H; =1 

X[l, if a PMU installed at bus i 
0, Otherwise 

f(X) is a vector function representing the constraints!is a 
vector whose entries are all equal to 1. 

PMU PLACEMENT RULES 

The following PMU placement rules are as follows: 

Rule 1: Assign one voltage measurement to a bus where a 
PMU has been placed, including one current measurement to 
each branch connected to the bus itself. 

Rule 2: Assign one voltage pseudo-measurement to each 
node reached by another equipped with a PMU. 

Rule 3: Assign one current pseudo-measurement to 

eachbranch connecting two buses where voltages are known. 
This allows interconnecting observed zones. 

Rule 4: Assign one current pseudo-measurement to 

eachbranch where current can be indirectly calculated by 
theKirchhoff current law. This rule applies when the 
currentbalance at one node is known, i.e. if the node has 
nopower injections (if N-l currents incident to the node 
areknown, the last current can be computed by difference). 

III. PMU FOR STATE ESTIMATION 

The PMUs receive signals from GPS satellites, and provide 
synchronized measurements from different locations to the 
desired destination, commonly known as the phasor data 
concentrator (PDC) [9]. The measurement data can be used 
for wide area monitoring; real time dynamics and stability 
monitoring; dynamic system ratings; and improvements in 
state estimation, protection, and control. 

A state estimator estimates the voltage magnitudes and phase 
angles at the buses by using the available measurements in the 
form of power injections, power flows, voltage magnitudes, 
or current through the branches. The voltage and current 
phasors measured by the PMUs can be used in a state 
estimator in two ways [10]: 

1) To increase the confidence in the available 
measurements.2) To replace the available conventional 
measurements. 

When a PMU is placed at a bus, it can measure the voltage 
phasor at that bus, as well as at the buses at the other end of all 
the incident lines, using the current phasor and the known line 
parameters. It is assumed that the PMU has a sufficient 
number of channels to measure the current phasors through 
the entire branches incident to the corresponding bus. 

It is to be noted here that the errors in the voltage and current 
measurements by the PMU and the transmission line 
parameters induce uncertainties in the estimated voltage 
phasorat the other end of the line [11]. The measurement 
uncertainty can be computed by the use of the classical 
propagation error theory, making use of the maximum 
measurement uncertainties provided by the manufacturer and 
the maximum measurement uncertainty in the line parameters 
[12]. The estimated quantity can be expressed as a function of 
the measurements and the parameters, X 1? ...,X n , i.e. 
y=f(xi,x 2 ...,x n ), where y is the estimated bus voltage 
magnitude or the estimated phase angle and are the 
measurements or parameters. The measurement uncertainty 
can be expressed as 
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Where Ax [is the maximum uncertainty associated with the 
quantity."*' r 

STATE ESTIMATION 

A phasor measurement unit (PMU) [7] equipped with a GPS 
receiver allows for synchronization of measurements, 
yielding accurately measured and time-stamped voltage phase 
angles.. The general conclusion is that PMUs [6] have greatly 
improved observability and accuracy of voltage angle 
estimates. Despite some opinions to the contrary, PMUs will 
not make state estimation obsolete even if they are available at 
every bus in the system. As we know, measurements are not 
perfect; thus a redundant set of measurements will still be 
needed in order to identify bad data. All of these 
measurements can be considered dynamic since snapshots are 
performed every few seconds. The status of the assets (line 
status, breaker status etc.) as well as network parameters can 
be considered as static measurements. The network topology 
processor in Fig. 2 determines the topology of the network 
from the telemeter status of circuit breakers. Having an 
observable set of measurements is a necessary, although not 
sufficient condition, for EMS computer applications. While it 
is desired, coordination across the network quite often does 
not happen in real-time. The reasons for not having real 
time-model are varied. While many control and monitoring 
functions is computer based, there are still functions handled 
by telephone calls between the system operator and utility 
control centers. It is a well known fact that control room 
technology is behind today’s state-of-the-art In the IT world 


System Measurements 



Figure 2: State Estimation block diagram 

HE [5,8] state estimator is one of the more important utility 
network management applications that could benefit from 
synchronized phasor measurements. Improvement in state 
accuracy could benefit from the more accurate and more 
frequent time-tagged information of voltages and currents 
from phasor measurement units (PMUs). The improved state 
brings forth improved accuracy and greater confidence on 
real time security margins. This allows operators to operate 
the system in more economical and secure fashion. For 
example, more confidence on calculated security margin 


would eliminate a lot of uncertainties associated with its 
value, thereby allowing larger power transfers, or provide a 
more rational basis for refusing transactions. Prior art in 
accounting for phasors in classical state estimation requires 
appending the measurement vector with bus voltage angle 
measurements. Classical state estimation (SE) with 
synchronized phasor measurements in a non-invasive fashion. 
The classical SE measurement vector, (the P and Q flows, and 
voltage magnitudes); rather they are integrated with the 
calculated SE estimates by means of a linear state estimation 
equation post priori. This proposed integration model 
improves on prior phasor-in-state estimation methods by 
accounting for the full set of current phasor data. This hybrid 
weighted least squares (WLS) and linear SE model could also 
be performed in wide area measurement system platforms and 
does not require communicating PMU data via SC AD A 
channels [4]. 

State Estimation [5, 8] is the process of assigning a value to an 
unknown system variables based on the measurements 
obtained from the system. State estimator has been widely 
used as an important tool for online monitoring, analysis and 
control of power systems. It is also exploited to filter 
redundant data, to eliminate incorrect measurements and to 
produce reliable state estimates. Entire power system 
measurements are obtained through RTU (Remote Terminal 
Unit) of SCADA systems [4] which have both analog and 
logic measurements. 

IV. LINEAR STATE ESTIMATOR 

PMU measurements may be included by a slightly different 
formulation of the traditional non-linear weighted least 
squares or they may be taken into consideration after a 
preliminary system state has already been determined. Even a 
small number of these precise measurements can weigh 
heavily on the accuracy of the overall state of the system. 
However, a true application of PMU technology to state 
estimation would have all of the traditional measurements of 
real and reactive power injections and current and voltage 
magnitudes replaced by bus voltage phasors and line current 
phasors. If only PMU measurements are used, there are also 
no complications from the use of both polar and rectangular 
values in the state estimation process, as would be done when 
including PMU measurements in traditional state estimators. 
If a state estimator could function with only PMU 
measurements as inputs then many issues associated with 
traditional state estimators could be resolved. Because PMUs 
are synchronized with GPS, the problem of scan time 
becomes irrelevant. One could imagine looking at the state of 
the power system with a traditional state estimator versus one 
which used explicitly synchronized PMU data and comparing 
it to putting on a pair of spectacles for the first time and finally 
having a focused view of the world. Once the problem of scan 
time has been erased, the only issue of time is the 
communication and computational delay between the 
collection of the measurements and the employment of useful 
information for decision making by the operation and control 
applications. Additionally, when using PMUs as metering 
devices, the state of the system is actually being directly 
measured. However, estimation is still necessary for including 
redundancy and bad data filtering. Because of this, the 
placement of the PMUs is critical for achieving a fully 
observable system with a sufficient amount of measurement 
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redundancy. The redundancy can be gained not only from the 
line currents (which are linearly related to the system state) 
but from redundant bus voltage measurements. 

In this section presented linear state estimation from a single 
phase perspective. First, the basic mathematical concept of 
linear state estimation will be presented on a pi-model 
transmission line to demonstrate how voltage and current 
measurements are related to the system state. 

Mathematical Modelling of 5-Bus System 

Simple 5-bus system used as an example network in this work 
is given below: 


Busl Bus 2 Bus 3 



Figure 3: Simple 5-bus network 


Br anch impedances of example network: 


From Bus 

To Bus 

R(pu.) 

X(p.u.) 

1 

2 

0.02 

0.06 

1 

4 

0.01 

0.08 

2 

3 

0.05 

0.10 

2 

5 

0.02 

0.07 

3 

5 

0.01 

0.07 

4 

5 

0.03 

0.09 


Table 1: Branch impedances of example network 
The bus admittance matrix will then take the following form. 
The rules for populating the admittance matrix of a network 
are derived from Kirchhoff’s laws of current injections into a 
node. Sparing the derivation of the equations, there are two 
simple rules that can be used to construct the network 
admittance matrix by inspection. 


1. The /7th element of the admittance matrix is the sum of the 


admittances of all of the lines connected to bus i. 


2. The ijth element of the admittance matrix is the negative of 
the admittance connecting bus i to bus j. 

The shunt impedances from each of the components are 
added to the diagonal elements corresponding to each of their 
respective busses in form of shunt admittance Yjj. 

Y 12 =l/z 12 (1) 


Zi 2 -Ri 2 + X 12 (2) 


Z 12 —0.02+ 0.06i 
Y 12 =l/Z 12 = -5.0+15.0i 

yio =gi + jbi (3) 

yj°=gj+jbj (4) 

Bus 1 is not directly connected with Bus 3. 
Then, 

y 13 = o 


For all other branches we can calculate the value of 
admittance 



V 



Y a ■ 

•• 


' V 1 

/ = 

i 2 

= 

*21 

y 22 

■- Y 2N 




- C iV- 



Y,Y2 

■■ Y j¥jY J 




the network admittance matrix Y, 


Y = 


■ 6.5 - 27,3i -5.0 +15.0/ 
-5.0 + 15.0i 12.7-36.2i 
0 -4.0 + 8.0i 

-1.5+12.3/ 0 


0 

-4.0 + 8.0i 
6 . 0 - 22 . 0 / 
0 


0 -3.7H3.2i -2.0fl4.0i 


-1.5 f 12.3i 0 

0 -3.7 f 13.2i 

0 -2.0fl4.0i 

4.8 - 22.3i -3.3 f lO.Oi 
-3.3 f lO.Oi 9.1 - 37.2/- 


V. RESULTS 

There are four main database files that represent the system 
model: the series impedance input file, the shunt susceptance 
input file, and the voltage measurement and current 
measurement tables. These four files contain all of the 
necessary indexing and topology information necessary to 
construct the system matrices discussed in the previous 
sections. Additionally, the formatting of the input files which 
contain the raw measurements from the field must be 
constructed in a predictable way such that the computer 
program will interpret them in the same way each time it is 
read. The output of each of the estimator applications must 
also be constructed in a similar manner so that the 
human-machine-interface (HMI) can present the information 
to the user in a valuable and accurate format. 

STATE ESTIMATOR TESTING RESULTS 

The desired function of any state estimator is to filter out the 
measurement error using redundancy and knowledge about 
the measurement errors. The testing script used to test the 
Matlab application of the three phase linear state estimator 
using the data created by the test to demonstrate this in several 
ways. 

1. The true value and the estimated value of a single state 

variable are plotted on the complex plane for each 
iteration of the testing procedure. The raw voltage 
measurements corresponding to that particular state 
variable are also included on the same plot. This 
type of plot shows the ability of the state estimator 
to filter out measurement error and one in on the 
true value of the state variable. 

2. A different plot which shows the estimators ability to 

filter out measurement error from the entire state 
vector. These plots only show results from a single 
iteration. Since the true value of the system state is 
known from the load flow, it is compared to the 
corresponding voltage measurements and estimated 
state vector and plotted to show the effects. This is 
done in four different plots: a plot showing the 
effects on the real part of the state variable, a plot 
showing the effects on the imaginary part of the 
state variable, a plot showing the effects on the 
magnitude of the state variable, and a plot showing 
the effects on the angle of the state variable. 

Figure 4 and Figure 5 are two plots showing a single state 
variable over 100 iterations of the testing script. This 
particular state variable corresponds to a substation which is 
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monitored by a PMU. Figure 4 shows in green the true value 
of the state variable and in blue the estimated values of the 
state variable for 100 iterations. Figure 5 shows in blue the 
estimated values of the state variable and in red the voltage 
measurements that correspond to this particular state variable. 
Both of these plots are from the same set of iterations so they 
can be visually compared to each other. These plots show 
very clearly the ability of the estimator to filter out 
measurement errors and determine the best estimate of the 
system state. It can be seen from Figure 6 that the true value of 
the state variable differs from the estimated value. Figure 7 
and Figure 8 are the same type of plot but show a state 
variable which corresponds to a substation which is not 
monitored by a PMU. 



Figure 4: State Estimator Output for 100 Iteration 



Figure 5: State Estimator Output and voltage measurements 
for 100 Iteration 



Figure 6: State Estimator Output for 100 Iteration (w/o 
PMU) 



Figure 7: State Estimator Output for 100 Iteration (w/o 
PMU) 

The next set of plots show the effects of the estimation 
process on the entire state vector for a single iteration of the 
testing script. The absolute value of the difference between 
the voltage measurements and the actual state vector are 
shown in red. The absolute value of the difference between 
the estimated state vector and the actual state vector are 
shown in blue. 



Figure 8: Comparison of Real Part of Output and 
Measurements 



Figure 9: Comparison of Imaginary Part of Output and 
Measurements 
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Figure 10: Comparison of Magnitude of Output and 
Measurements 



Figure 11: Comparison of Angle of Output and 
Measurements 

It can be seen from these plots that the state estimator can 
effectively filter out measurement normally distributed 
random additive errors and estimate with a good degree of 
precision the system state. In the case of all four comparisons, 
the estimator has the ability to come within approximately 
0.002 of the true value of each of the parameters. 

The database files used by the topology processor application 
include a current measurement lookup table, a breaker 
location lookup table, and a bus/branch model describing the 
full network. This section describes the composition of each 
of these database files and provides examples. The current 
measurement lookup table for the topology processor is not as 
detailed as the current measurement lookup table for the state 
estimator. This lookup table only contains the current 
measurement numbers and the location of each of those 
current measurements in the network. The location is given as 
the line number the measurement originated from. Below is an 
example of the current measurement lookup table 2 for the 
topology processor application. This table corresponds to the 
simple 5 bus system in Figure 3 


Current measurement location:- 
Lookup table = [ 


Current 

Mess 

Line 

Number 

1 

1; 

2 

1; 

3 

2; 

4 

2; 

5 

3; 

6 

3; 

7 

4; 

8 

4; 

9 

5; 

10 

5; 

11 

6; 

12 

6]; 


Table 2: current measurement location 

The breaker location lookup table is slightly more 
complicated. It is organized by the line numbers of each 
branch in the network. Following each line number is a series 
of breaker numbers which all must be “open” in order for the 
line to be out of service. For this application it has been 
assumed that no more than four breakers would be required to 
open any particular line. Other breaker configurations in the 
substation could yield lesser numbers. In this case, the 
unfilled slots in the input file should be filled with zeros. 
Below is an example of the breaker lookup table for the 
topology processor application. This also corresponds to the 
simple 5 bus system in Figure 3 
Breaker Lookup Table:- 

Loo kup table = [ ____ 


Line Number 

BS1 

BS2 

BS3 

BS4 

1 

1 

2 

4 

5; 

2 

4 

6 

7 

8; 

3 

2 

3 

10 

ll; 

4 

5 

6 

13 

14; 

5 

8 

9 

13 

15; 

6 

11 

12 

14 

15]; 


Table 3: Breaker Lookup Table 


Finally, it is valuable for the topology processor to have 
access to the bus/branch model of the full network. This is 
contained in another database file and is formatted as follows. 
Full Network Bus/Branch Model:- 
Full network = [ _ 


Line 

Number 

From 

Bus 

To 

Bus 

Status 

1 

1 

2 

l; 

2 

2 

3 

l; 

3 

1 

4 

l; 

4 

2 

5 

l; 

5 

3 

5 

l; 

6 

4 

5 

l]; 


Table 4: Full Network Bus/Branch Model 
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Conclusion 

The ability of the state estimator to achieve a high level of 
efficiency and numerical robustness is of paramount 
importance in today’s electric utility industry. A robust 
algorithm must be globally convergent(convergent from any 
starting point), and able to solve in practice both 
well-conditioned and ill-conditioned problems. 

The objective is to provide a more reliable and robust state 
estimator, which can successfully cope with all kinds of errors 
(bad data, topological, parameter) faced in power system 
models. 

An important advantage of this method is non-iterative which 
reduces the computation time and suitable to implement in the 
near future. In view of extension in this work, analysis can be 
done taking into account the loss of measurements, failure of 
PMU and effectiveness of PMU. With the increasing use of 
computers and digital information in power systems, the 
availability of information about a system has become easier. 
Online data, artificial neural networks, fast computations and 
real time control are beginning to take over the theoretical and 
analytical solutions. However, the conventional methods of 
power system analysis still provide the element of theory to 
the latest methods of computations and estimations. Within 
this environment, there will be an increasing need of more and 
more such application oriented studies which can help 
minimize the resources and computation time. This study was 
intended to show that a small amount of real time information 
with limited resources can be expanded into a bigger picture 
of the whole system with simple relationships as derived 
herein. 

In three phase linear state estimator feedback of data comes 
from PMU with help of communication system. When there 
is any problem find found in PMU or communication system 
it then there is no feedback of data. Thus there should be one 
more source of data is required. At present we get non-linear 
feedback of data with help of SCADA. By using these two 
methods, in future we can develop a logical state estimator for 
overcome these type of problem 
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